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Abstract 



Dark soliton excitations are shown to exist in single wall carbon nan- 
otubes (SWCNTs). At first, the nonlinear effective interatomic potential and 
the difference equation for longitudinal lattice displacement are obtained for 
the SWCNTs by expanding Brenner's many-body potential in a Taylor series 
up to fourth-order terms. Then using a multi-scale method, for short wave- 
length lattice excitations the equation of motion of lattice is reduced to the 
cubic nonlinear Schrodinger equation. Finally, the dark soliton solutions and 
relevant excitations in the SWCNTs with subsonic velocity are discussed. 
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Carbon nanotube (CNT) is the name of ultrathin carbon fiber with nanometer-size di- 
ameter and micrometer-size length and was accidentally discovered by Sumio Iijima ||T| in 
1991. As a novel and potential carbon material, CNTs have received a great deal of attention 
and their many unique properties, such as structure [|7-TI[], electronic property fl2|-|14 



superconductivity [0-0 and elementary excitations [|T8|-|2"2]] , have been investigated. The 
structure of the CNT consists of an enrolled graphite sheet, and can be classified into either 
multi-wall or single wall CNT (MWCNT or SWCNT) depending on its preparation method. 
The smallest SWCNT with a diameter of 4 A was discovered in last year fl23| , |24" . 



In recent years, there has been considerable efforts for studying the nonlinear localized 
excitations in lattices [p5|-p7|. The nonlinear localized excitations can transfer energy and 



be involved in various processes of interest. There have been some speculations on the role 
played by solitary excitations in heat transfer, polymer destruction, and other processes oc- 
curring in molecular systems. The dynamical properties of the CNTs are of great interest due 
to their potential for useful practical applications |28j . The longitudinal soliton excitations 
governed by Korteweg-de Vries (KdV) equation in the CNTs using Brenner's many-body 
potential have been investigated in Ref. P2| . Such solitons are obtained under the assump- 
tion of a long wavelength approximation. It is known that, in one-dimensional (ID) lattices, 
for a weekly nonlinear lattice excitation with a large spatial extension, its amplitude (or 
envelope) is governed by the nonlinear Schrodinger (NLS) equation when the excitation is 
a short wave wavepacket p9|| . The nonlinear excitations in this case are envelope solitons. 
It is interesting to consider the envelope soliton excitations in the SWCNTs. 

In this work we investigate the nonlinear effects in SWCNTs giving rise to dark envelope 
solitons using Brenner's empirical many-body potential for carbon systems p0| . At first, 
we introduce a quasi-one dimensional lattice model for armchair (tube (m,m) ) CNTs in 
an anharmonic approximation of the simple analytical Brenner's potential. Then we derive 
a NLS equation by using a multiple-scale approach. Finally, we discuss the dark soliton 
solutions of the NLS equation and the physical relevance of these excitations in the SWCNTs. 

The SWCNTs are highly anisotropic and ordered objects. The surface of the SWCNT 



is formed by a graphite sheet folded into a cylinder with bond lengths and angles differing 
slightly from graphite on account of the strain induced by the folding. The interaction 
between adjacent atoms i and j in the SWCNTs can be described by the Brenner's many- 
body potential j22"fl , which reads 

4 \ - BijV£, (i) 

where V$ and are, respectively, exponential repulsive and attractive terms: V$ = 
27.27 exp[-3.28(ry - 1.39)] and V$ = 33.27 exp[-2.69(ry - 1.39)]. By represents an en- 
vironment dependent many body coupling between atoms i and j containing geometric 
information associated with the system. The total energy of the system is obtained by 
summation of Eq. ([J) over all bonds. The energy and distances are measured in eV and 
angtroms, respectively. 

We take the cylindrical coordinate system (R, $, z) and align the tube axis as the z 
axis so that 2m atoms in the SWCNT layer have the same z coordinates. The equilibrium 
distance between layers, l = 1.26A for a (5,5) tube. 

Let us consider a cylindrically symmetrical disturbance in the SWCNT geometry as done 
in Ref. [Q, in which all 2m atoms in the nth layer have identical z and radical displacements 
from their equilibrium positions Z® and R°: Z n = Z® + ( n and R n = R° + p n , where Z n 
and R n are the perturbed coordinates and ( n and p n are the displacements from equilibrium 
positions. Then the coordinates of ith atom in nth layer are (-R + p n , Z® + ^ n ) for weak 
nonlinear excitations the ( and p are much smaller than the characteristic length scale. So 
the interatomic potential can be expanded in a Taylor series and the terms up to fourth 
order retained: 

E = E + J2 2™E n (2) 

71=1 

where the E is the ground state energy of relaxed SWCNT, E n = \{E\ + E 2 + E 3 ) is the 
energy of an atom in the nth layer due to atomic displacements, where E\, E 2 and E 3 are 
the energies of bonds emerging from any atom in the nth layer. 
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We rewrite the Brenner's potential (]T|) by expanding it in Taylor series up to the fourth- 
order terms as 



E(r) = — (cr 2 exp[-cri(r - r )] - ai exp[-cr 2 (r - r )]) 

= {( r - r o) 2 - ^-a x {r - r ) 3 + /^O - r ) 4 , (3) 
2 3r 4rp 



where A = 27.27 eV, o x = 3.28A- 1 , a 2 = 2.69A" 1 and r = 1.39A and / = ^Afa - 
a 2 ) = 52.77eV/A 2 , a x = f (oi + a 2 ) = 4.15 and a 2 = f (of + <7i«7 2 + of) = 8.64. Thus 
Ei = E(n),i= 1,2,3 with 



{(R° + p„) 2 (cos $ m - cos $i) 2 + + p„) 2 (sin $ m - sin $,) 2 } 1/2 = r (l + pji? ), 



r 2 = K + (p n+1 - Pn f + (r 2 - Z 2 )(^^ + 80*) + 2 / o(Cn+1 - C n) + (Cn + i - Cn) 2 } 1/2 

and r 3 = r 2 (n + 1 — > n — 1). Expanding r 2 and retaining the terms to the third-order, we 
get 

r 2 = r + X 1 +X 2 + X 3 , 



where 



r 2 — I 2 l 
Xl = 2R°r ^ Pn+1 + Pn ^ + r^ n+1 ~ 

1 r 2 - I 2 (r 2 - I 2 ) 2 

X 2 = — (Pn+i - Pn) 2 + ^o^/WiPn - grWV + P " )2 

(/ " " ^(Pn+l + Pn)(Cn + l - Cn) + ^(1 - ^XCn+l " Q\ (4) 



2r 3 R ^n + x ■ ^ASnt. sn, ■ ^ 
X3 = ~~^3 ( (Pn+1 - Pn) + Jjjpy Pn+lPn + (Cn+1 ~ Cn) 



2R° 



(p n+ l + p n ) + Z (Cn+l - Cn) 



Then we have 



" ~ 4(fl») 2P " 6(fi»)3" lP " 8(fi») 1<l2P " 



i, x! + fXlX2 _ |^ + n \xl + xa, - Ixlx, + %x }) 



4 



If we neglect the radial degrees of freedom and reduce £ to the dimensionless unit ( n — > 

C„// , we get 

E n = ^(Cn+l - Cn) 2 ~ y (Cn+1 ~ Cn) 3 + ^"(Cn+1 ~ Cn) 4 , (6) 

where b = ^l 2 J = 68.84eV,p = -f+(§+ai)f = 3.14, q = I-(2a 1 +3)|+(i+2a 1 +a 2 ) J = 
2.99. The equation of motion for longitudinal displacements reads 

^ = [C„ +1 - 2( n + Cn-l] " P [(U " Cn) 2 " (Cn " Cn-l) 2 ] 

+g[(Cn + l-Cn) 3 -(Cn-Cn-l) 3 ], (7) 



where we have made the dimensionless transformation for time t as t' — Jb/ '(Mlfyt, M is 
the mass of carbon atom. 



Now we investigate how the narrow band wave packets evolve by nonlinear effect p9| , pT 
We look for a solution expanded in terms of a small but finite parameter e denoting the 
relative amplitude of the excitations: 



OC 



Cn = E e^tiJ (8) 

and each «M i s also expanded in terms of harmonics exp(U9 n ) with "fast" variable 9 n = 
knlo — out' as 

oo 

«J° = E u { :}(T,t n )e^. (9) 

l=— oo 

where r and £ n are slowly vary variables defined by r = e 2 t' and £ n = e(nlo—\t'), respectively. 
A denotes the group velocity A = ^ which will be determined later. It is noted that the 
reality conditions 

<] = (io) 

should be satisfied for all z/ and I. 

Substituting Eqs. (§) and (|9[) into Eq. (0) and equating the coefficients of various powers 
of e to zero, we get 



E{-/V + 4sin 2 ^} M i°)e^ = 0, (11) 
l / 



i \ 



^ a iT,7'W„„( )„(°)^(W)*» 



«n,z e 



= PEE 8 « sin 2 ^ sin ZU<!«S e ' (WOfl " > ( 12 ) 
z Z' A 

d 2 u (0) g u (0) du (l) 

{(A 2 — I 2 , coslkln) - 2ilu>— + 2i(luX — L sin /fc/ ) - "' 

+ (_/V + 4W^)4i}e^ 



= pEEI 8 * sin2 sin ^ fe ^o(wSl u Si' + 

Z Z' z 

+8 / sin 2 ^ cos I'klou^J -^f + 4/ sin lkl sin J' W ^«S }e l « +l '^ 
-9EEE 4sin2 — {(e u ' kl ° - l)(e«" Wo - 1) + (e a ' feZ ° - 1)(1 - e - fl " Wo ) 

Z Z' Z" 2 

+ (1 - e"^)(l - c-^)}!^^ uS^'^- (13) 

and the equations corresponding to higher order of e. To derive these equations, we have 
expanded the terms u„±i \ in Taylor series of eIq, because £ n ±i = £ n ± eIq. 

From above equations, we obtain the dispersion relation u 2 = 4 sin 2 ™, here k and a; 
being respectively the wave number and the frequency of a lattice wave, and the following 
equations for the envelopes of the lattice wave 

vQ = 0, for |/| > 2 (14) 
d 2 u { nl _8p d (o) 



■ att n,l f ° U n,l . r<J 2 , 2/,i 2\l I (0) |2 (0) j OU n,0 (0) / 1r n, 

»-^T = | w ^f + <*>{-gw + p (4 - w )}K/ 1 | </ x - pZow-^-<,i. (16) 
The Eq. fll5|) is easily integrated once and written as 

dUnfi _ 8 P (\ (0),2 



^{l«S 2 -C}, (17) 
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where the — is an integration constant. Therefore from Eq. fll6|) we obtain the single 
equation to determine 

,OU nl L Q O U nl d 2 It a , 2\ 1 I (0) 1 2 (0) i o 2 /"r (0) / 10 \ 

*-^T = ^"T^ - +a; \2 qUJ ~P ( A + UJ )j Kil U n,l + ^<i- (18) 

By replacements of variables, we obtain the standard cubic nonlinear Schrodinger equation 

^ ld 2 d) , a 



<9z 2 <9x 2 



(19) 



where = exp(i8p 2 uj J Cdr)u^{, z = ut, x = f-£ n and 

« = (P 2 - ^)^ 2 + V- (20) 

Here k is always positive, because p 2 — |g = 5.37 > 0. It is known that the reversal of the 
sign of k leads not only to a change in the physics picture of the phenomena described by 
Eq. flT9|), but also to a considerable restructuring of the mathematical formalism necessary 
for its solution. 

Now we investigate the modulational stability and solitons |!2]|!3|. The NLS Eq. ( |T9"D 
has the simplest solution in the form of a continuous wave given by the expression, 

= oe <"*-H/fc = (I« 2 + (21) 

where <fto is a constant. Let us investigate the linear stability of the exact solution, Eq. (plD , 
against small perturbation as 

<f>=((f> + x)e i{ax+Pz+, "\ (22) 

where the function \ an d derivative of the phase if) are assumed to be small. Substituting 
Eq. fl2"2|) into Eq. (|T9|) , we look for the solution to these functions in the form, x, ?/> ~ 
exp(iQz — iQx), then obtain the dispersion relation as 

Q 2 (K(f)l + ^Q 2 ) = (aQ + n) 2 (23) 

which means that the small excitation are always stable when k > and the small amplitude 
waves can propagate along the background. As a result, Eq. (|19"D has soliton solutions in 



the form of localized "dark" pulses created on the continuous wave background. The NLS 
equation with the boundary conditions \<f>\ — > (fi , at x — > ±00, is exactly integrable by the 
inverse scattering method 



(cos p tanh O + i sin ip)e 



iax+ifiz 



(24) 



with = t/k.(Po cos <p[x — xo + (a + a/k0o sin tp)z\ , where </>q, ip, a and xq are free parameters. 
The 0o is the amplitude, and p corresponds to the total phase shift across the dark soliton, 
2<p. Then we obtain 



Ap 



(o) 

U n.O 



a/k^q COS if 



— C)6 — o cos yjtanhB 



(25) 



If we take C = cfl for simplicity, then we obtain dark soliton excitations for the longitudinal 
displacement ( n up to the lowest order as 



( n = A tanh 9 [cos X 



2p 



K 



Aq tan p sin X 



(26) 



with 



X = {k + —a')nlo 



nAl 



2a cos h (2 — a 

2 2 COS 2 <y2 COS 2 <£> 



+ — — ) sm — 



M/ 2 ' 
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K r fcZ 

t 2 



a 



— ylotanv?) sin — 



where A n 



io£ cos y> is the rescaling amplitude and a' = ea is the rescaling parameter. 



The velocity of the dark soliton is given by v g = [cos ^ — (a' + tamp) sin ^]v sound , 



where v soun d = yb/M = 22km/sec is the longitudinal sound velocity, which means the the 
velocity of dark soliton in the SWCNT is subsonic. 

To summarize, We have obtained the nonlinear effective interatomic potential for the 
SWCNTs by expanding Brenner's many-body potential in a Taylor series up to the fourth- 
order terms and the difference equation for longitudinal lattice displacements. Using the 
multi-scale method, for short wavelength excitations the equation of motion of lattice has 



S 



been reduced to the cubic nonlinear Schrodinger equation. We have discussed also the 
envelope dark soliton excitations in the SWCNTS. The dark soliton is the excitation of a dip 
excited on the continuous background of a high frequency harmonic oscillation of longitudinal 
displacement, and its velocity is subsonic. We note that the KdV soliton obtained in Ref. 
is a weak nonlinear excitation valid only in a long wavelength approximation with a 



supersonic propagating velocity. 
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